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ABSTRACT 

Aims. To probe cosmological fields beyond the Gaussian level, three-point statistics can be used, all of which are related to the 
bispectrum. Hence, measurements of CMB anisotropics, galaxy clustering, and weak gravitational lensing alike have to rely upon an 
accurate theoretical background concerning the bispectrum and its noise properties. If only small portions of the sky are considered, 
it is often desirable to perform the analysis in the flat-sky limit. We aim at a formal, detailed derivation of the bispectrum covariance 
in the flat-sky approximation, focusing on a pure two-dimensional Fourier-plane approach. 

Methods. We define an unbiased estimator of the bispectrum, which takes the average over the overlap of annuli in Fourier space, and 
compute its full covariance. The outcome of our formalism is compared to the flat-sky spherical harmonic approximation in terms of 
the covariance, the behavior under parity transformations, and the information content. We introduce a geometrical interpretation of 
the averaging process in the estimator, thus providing an intuitive understanding. 

Results. Contrary to foregoing work, we find a difference by a factor of two between the covariances of the Fourier-plane and 
the spherical harmonic approach. We argue that this discrepancy can be explained by the differing behavior with respect to parity. 
However, in an exemplary analysis it is demonstrated that the Fisher information of both formalisms agrees to high accuracy. Via 
the geometrical interpretation we are able to link the normalization in the bispectrum estimator to the area enclosed by the triangle 
configuration at consideration as well as to the Wigner symbol, which leads to convenient approximation formulae for the covariances 
of both approaches. 
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1. Introduction 
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As the concordance model of cosmology becomes more and more consolidated, the focus increasingly turns towards probing ef- 
fects beyond the standard paradigm, such as non-Gaussian initial conditions or the evolution of the large-scale structure in the 
highly non-linear regime. To lowest order, these effects can be measured by three-point statistics of the underlying fields, all of 
] which are related to the bispectrum. Hence, work in both t heory and observations concerning the bispectrum and its noise proper- 
i ties has been undertaken for CMB measureme nts (e .g. lHull2000l : ICoorav et alJl2008h . galaxy clustering surveys (e.g. IScoccimarrol 
' 2000t IScoccimarro et alJ 2001 Sefusatti et alJ 2006t). or. more recently, weak gravitational lensing on cosmological scales (e.g. 
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^ ! iBernardeau et al.ll2002bHJarvis et al.ll2004tlTakada & Jainll2004l) . 
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While theoretical computations at the bispectrum level are already considerably more demanding than for second-order statistics, 
| this does apply even more so to the bispectrum covariance, which is a six-point statistic. On the full sky calculations are done by 
• • . expanding the signal into spherical harmonics. If only small angular scales are considered, it is often more convenient to use a 
flat-sky approximation and work in terms of Fourier amplitudes. In the case of weak lensing the flat-sky limit is appropriate for 
practically all applications because signal correlations can only be measured up to separations of a fe w degrees. 

Although other approaches exist in the li teratu re (e.g. lMatarrese et al.ll 19971: ISefusatti et al1 l20061. a lot of work is done within 
d ' a flat-sky spherical harmonic formalism (Hu 2000), which suffers - at least formally - from drawbacks. For instance, the resulting 
flat-sky expressions are valid only for integer arguments and thus for a bin width of unity, whereas it is desirable to evaluate the 
bispectrum and its covariance at real-valued angular frequencies and e.g. a logarithmic binning. The formulae still contain Wigner 
symbols whose physical meaning within a flat-sky consideration remain obscure. As the spherical harmonic expansion can only 
be done on the full unit sphere, the finite size of the survey at consideration is usually accounted for by multiplying a factor, 
containing the sky coverag e, by hand . Moreover, the accuracy of some of the approximations in the transition between full sky and 
two-dimensional plane (see Hu 2000) is uncertain. 

This work aims at clarifying the derivation of bispectrum covariances in the flat-sky limit. We attempt to do so by presenting a 
detailed calculation which is purely based on the two-dimensional Fourier formalism, followed by a comparison of this approach 
with the flat-sky spherical harmonic results in terms of their covariance, the behavior under parity transformations, and the informa- 
tion content. Moreover, we provide further insight and illustration by establishing relations between Wigner symbols, the averaging 
process in the bispectrum estimator, and a geometrical view. 
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The outline of this note is as follows: In Sect.[2]a bispectrum estimator is defined and shown to be unbiased. Section0introduces a 
geometrical interpretation, which is then applied to deal with the issue of degenerate triangle configurations. In Sect.|4]the covariance 
of the estimator defined beforehand is computed. The result is compared with the spherical harmonics approach and demonstrated 
to be equivalent in terms of information content in Sect.0 To explain the differences between the covariances, we also discuss the 
treatment of parity in both formalisms. We summarize our findings and conclude in Sect. [6] To avoid confusion, we r efrain from 
using the term 'flat sky' in the following, but refer to our formalism as 'Fourier-plane' and to the approach as e.g. given in Eu|{2000) 
as 'spherical harmonic' (both are flat-sky approximations). 



2. Bispectrum estimator 

We consider a continuous, two-dimensional random field g with mean zero, which is characterized by its complex Fourier amplitudes 
g(f), where I denotes the angular frequency vector. Throughout, it will be assumed that this field is statistically homogeneous, i.e. 
invariant under translations, and statistically isotropic, i.e. invariant under rotations. In a cosmological context g could for instance 
represent the temperature fluctuations of the CMB, the num ber density contrast of galaxy surveys, or the weak lensing convergence. 

In what follows we will largely follow the approach of Joachimi et al. (2008), assuming likewise measurements in a compact, 
contiguous survey of size A. We will restrict our considerations to an angular extent much smaller than the size of the survey, i.e. to 
I » 7r/9 max , wh ere fl max is the maximu m separation allowed by the survey geometry. Boundary effects due to the finite field size, as 
e.g. discussed in Joachi mi et alJ {2008) for the second-order level, can then be safely neglected. 

Furthermore, we will not explicitly consider additional noise terms due to the discrete sampling of the continuous field g, for 
ease of notation. To account for these shot noise or, in the case of weak lensing, shape noise terms in the covariance, they c an simply 
be added to the second-order measures, so in this Fourier space approach, to the power spectra (e.g. Kaiser 1998k IHtj|I 1999h . 

Note that the galaxy ellipticity, and not the convergence k, is the direct obser vable in weak lensing. However, in absence of shape 
noise and for ( » 1, the estimators in terms of the galaxy ellipticity, as given in Joachi mi et alJ d2008l) . can be re-written directly in 
terms of k. Thus, without loss of generality, one can consider the convergence as the observable that the estimator is based on. 

For a statistically homogeneous and isotropic random field one defines the bispectrum as 

igitx) g(l 2 ) gih)) = (2nf Sff(ii + l 2 + h) B(AJi, h) , (1) 

where 6^\f) is the two-dimensional Dirac delta-distribution. It ensures in ([TJ that the three angular frequency vectors form a 
triangle. For the assumed properties of g the bispectrum has three independent components, for which we have chosen the triangle 
side lengths \((\ = £j. For the absolute values of ( the triangle condition translates into the requirement \£i - £ 2 \ < £3 < £\ + £ 2 or 
equivalently for any permutation of th e 

Similarly to I Joachimi et alJ d2008l) . we construct an estimator of the bispectrum by averaging configurations over annuli, where 
here one has the complication of allowing only those combinations of angular frequency vectors that form a triangle. The area of an 
annulus with mean radius l{ is given by 

A R (ft) = InlMi (2) 
with the bin size A€,. Then we define the estimator 

Kh,~hM := ^A" 1 (*ub,h) f 4lhA Tlk\ T^^D 2) (A+f2+ ^3) ^1)^2)^3), 0) 
A v ; Ja^) A r (({) J At <to A R (( 2 ) Ja 8 (« A r (€ 3 ) 

where A is a function that is related to the fraction of angular frequency combinations allowed by the triangle condition. It is 
defined such that (0 is unbiased, its explicit form being calculated below. Note that this bispectrum estimator is invariant under any 
permutation of its arguments since A is symmetric as will be shown below. 

In the following, we demonstrate that (0 is unbiased by computing the ensemble average, 

m,k.ty = &£A- 1 (i u k.h) f 4ik f 4ik f Pr(2n) 2 {^({ l+ e 2 + { } )) 2 B(£ u £ 2 j } ) (4) 

N ' A > Ja,^) A r {£ x ) J Al! Cc 2 ) A r (£ 2 ) J Ag( « A R (h) VD ' 

= (2^a- 1 (a,I 2) I 3 ) f -f^-f ^-^(e l+ e 2 + e,)B( (l ,£ 2 j 3 ). 

JAsih) A R\£\) JA R (t 2 ) A R (£ 2 ) jAsih) A R (£ 3 ) 

In the first step the definition of the bispectrum (0 w as inserted. The appearance of a squared delta-distribution requires taking into 
account the finite survey size. As shown in Joach imi et al.l (120081) . one can identify 

which results in the second equality of (0|. 

Since the bispectrum only depends on the magnitudes of the angular frequency vectors we can perform the integrations over the 
polar angles of the ^-integrals. If tp(. denotes the polar angle of f,, one gets 

r*2jl r*2n r*2jr r*2n r*2n r*2jz r* f\^Q 

dip tl &p i2 d^6% > (e l +e 2 + e } )= d^, d^ 2 d^ 3 e 1(f ' + ^ )fl (6) 

Jo Jo Jo Jo Jo Jo J \ Z7: ) 

= J ^_(2nfj ((y8) MC 2 0) MhO) = (2n) 2 J AG 6 Wrf) J {€ 2 ff) Wrf) = 2jt A(£ u £ 2 ,£,) . 
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Fig. 1. Sketch of the annuli and their 
overlap for fixed t \ . The region of over- 
lap is approximated by the shaded par- 
allelograms. Note that due to mirror 
symmetry a second shaded area, related 
to the triangle t\, £L C v contributes as 
well. 



After inserting one possible representation of the delta-distribution in the first equality, we have made use of the definition of the 
Bessel function of the first kind of order 0, 



f 
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The result of the integral over three Bessel functions is taken from Gradsht evn et alJ d2000l) . formula no. 6.578.9, where we have 
defined 

{ else 

i.e. if ^3 are chosen such that they can form the sides of a triangle, then A -1 is the area of this triangle. Hence, © represents 
th e de f ining e quation for A. The set of integrations (O is also performed within the spherical harmonic approach, see the appendix 
of |Hu| ((2000), with a different result, which will be investigated in Sect. 15. II Note furthermore that A (l\,t2, ^3) = in case the 
angular frequency vectors are collinear or equivalently, if t t + ij = 4 for some combination i,j,k e {1,2,3}. At the same time, 
the bispectrum is non-zero for these degenerate triangle configurations, see (Q]i. For the time being, we exclude degenerate triangles 
from the derivation, but develop a treatment for these cases in Sect. 13.21 
Inserting (0 into one obtains 

{Ji+l/l&e p ^ 2 +l/2A£ j/ p w\ + l/2Af j» f 

(B(e u t 2 M) = (2nfh- l ((uh,h) A(A,4,^)B(A^2,«. (9) 

X ' X ' J^-VlKt A R ((i) J?,-l/2Af A R (( 2 ) J? 3 -l/2Af A R (i 3 ) 

Analogous to the derivation at the level of second-order statistics (I Joachimi et al.ll2008h we assume now that the annuli are thin 
enough such that A within the integral, evaluated at the average ^-values, can be taken out of the integration. Applying in addition 
©, one arrives at 

(B(e u h,h))~ B(£iJ 2 ,{ 3 ) = B(£ u e 2 Ji) , (10) 



where in the last step the bin-averaged bispectrum was defined. Hence, OJ defines an unbiased estimator of the bispectrum. 
Following the restrictions on ©, this estimator is non-zero if the condition \l\ - l 2 \ < Ij, < (\ + £ 2 , or likewise for all permu- 
tations, holds. 



3. Averaging over triangles 

A central step in the construction of the bispectrum estimator (0 is the correct treatment of the averaging over annuli, given the 
triangle condition. This section provides an illustrative, geometrical interpretation of the averaging process and applies this view to 
a practical treatment of degenerate triangle configurations. 



3.1. Geometrical interpretation 

Without loss of generality consider t\ to be fixed. Due to the assumed statistical isotropy of the underlying random field the angular 
integration over tp^ is expected to simply reduce to an average over all directions of t\. Then the geometric situation in the Fourier 
plane can be seen as in Fig.Q] For a given triangle, composed of the mean vectors t\,t 2 , l 3 with lengths {, for i = {1, 2, 3), the annuli 
for ( 2 and {3 are shown. Due to the triangle condition, the average is not taken over the whole area of the annuli, but merely over 
the region that the annuli have in common. This area of overlap is well approximated by a parallelogram of size A\\ = A£ 2 Aij,/ sin or, 
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Fig. 2. Comparison of expressions for the overlap area of annuli. Left panel: Relative deviation of ([T2l > from the overlap area of the 
annuli as a function of angular frequency. The bin width is kept constant at A = 0.05. The solid curve shows results for l\,€2 = 200 
and varying (3, while the dashed curve corresponds to ~t\ = 200, £3 = 400 and varying i 2 . Center panel: Same as above for the case 
t\ = 200, {3 = 400 and varying I2, but with A = 1. The dotted line illustrates the deviation of ( fT2l . the dashed line the deviation of 
(TT~3T >. Right panel: Area of overlap for the case t\ = 200, {3 - 400 and varying 1%, with A = 1. The solid curve corresponds to the 
actual area, the dotted curve to (1121 1. and the dashed curve to (fl3T l, Note that due to A = 1 only the values at integer values of are 
relevant for the covariance calculation. 



where a is the internal angle of the triangle opposite t\. This relation can readily be computed from the geometry of the sketch and 
by noting sin a = sin(7T - a). 

The configuration is mirror-symmetric with respect to an axis through l\. Correspondingly, another area of overlap of the same 
size, which is connected to the triangle l\, C 2 , £' 3 , contributes as well. Noting that axis reflection is in two dimensions equivalent to 
the parity transformation, the averaging is performed over triangles of both parities. A detailed discussion on this issue is given in 
Sect.E2] 

As the angle a can also be related to the size of the triangle at consideration, A -1 = (1/2) £2^3 sin a, one finds the following 
correspondence of expressions: 

Jo 2n Jo 2n J 2n D V V A R ({ 2 ) A R (h) 

where the first equality is an immediate consequence of ©. To arrive at the last expression, we used (f2j). Hence, the angular 
integration over the delta-distribution yields the ratio of the area of overlap A over i ap , approximated by 2 Ay, and the product of the 
area of the annuli the ^-integrations (excluding the fixed t \ ) run over. This ratio is in turn proportional to the inverse of the area of 
the triangle spanned by the angular frequency vectors. Therefore, by placing a prefactor of A -1 in the estimator (0, one replaces the 
normalization by the area of the annuli with the effective area, over which the average is actually performed. 

Two approximations are involved in this picture. First, the shaded regions in Fig. Q] are approximated as parallelograms, which 
is a good assumption if the angle, at which the two annuli intersect, does not become too small. Moreover, the narrower the 
annuli, the less discrepancy between the area of the parallelogram and the actual overlap is expected. If the triangle approaches 
the degenerate case, where £2 and {3 eventually come to lie on l\, the area of overlap attains a more complex shape. In particular, 
the correspondence to the area of the triangle, whose inverse is divergent, does not hold anymore. Second, reconsidering (O, we 
have replaced the average of A over triangle side lengths by A, evaluated at the average side lengths. This approximation similarly 
breaks down for thick annuli and configurations in which a small change in the length of an angular frequency vector causes a strong 
change in the size of the overlap region, as is the case near degeneracy. 

In Fig.|2l we have plotted the relative deviation of 

A ov d erl a p = ^ ^2 A£ 3 A (l U i 2 , h) (12) 

from the actual area of the overlap region, which we calculated numerically. For simplicity, we assume a constant bin width 
A^2 = A^3 = A for all computations related to Fig. [2] For a small bin width A = 0.05, given integer steps in I, we find for the 
two configurations considered in the top panel that the approximation of the overlap area by parallelograms is excellent for the 
vast majority of triangle configurations. However, as expected, the deviation rises sharply when approaching the degenerate case. 
Changing to A = 1, i.e. the maximum meaningful bin width in this setup, the relative deviation is larger, but still very small except 
for triangles close to degeneracy. 
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Fig. 3. Sketch of the region averaged 
over in case of a degenerate triangle, 
again for fixed t\. The depicted trian- 
gle has side lengths l\, £2 + A£ 2 /2, and 
£3 + At?,/ 2. The shaded parallelogram 
approximates the region of overlap, the 
mirror-symmetric counterpart not be- 
ing shown. 



3.2. Degenerate triangles 

As discussed in the foregoing section, the approximations made in the course of the construction of the bispectrum estimator break 
down for degenerate triangle configurations. Equation ( fTTT i becomes invalid, the inverse area of the triangle A diverging. Yet, to 
be of practical use, it is necessary to extend the validity of (0 to the case of degenerate triangles. We do so by making use of the 
geometrical interpretation of the averaging process. 

Still keeping t\ fixed, consider the situation of a degenerate triangle as sketched in Fig. [3] Here, £\ -£2 + (3, while the depicted 
triangle has side lengths l\, (2 + A£ 2 /2, and £3 + A^/2. Again, we identify a parallelogram that serves as an approximation for the 
overlap of the annuli, although, as the sketch suggests, with considerably lower accuracy. The relation between the internal angle 
a of the triangle to the internal angle of the parallelogram n - a holds as before, so that one can derive an analogous formula to 
(TT2l . but with modified triangle side lengths. Symmetrizing this argument for all three angular frequency vectors, we propose the 
following formula to compute the area of overlap in the degenerate case: 



As is evident from Fig. [2 center panel, the relative deviation of ( fT3l from the true overlap area is still fairly small, but - 
unsurprisingly - noticeably stronger than for dl2l . The right-hand panel gives the size of the overlap area for values of £2 close 
to 200, which is the degenerate case. Note that since this plot was determined for A = 1, the values relevant for the covariance 
calculation are only those at integer I. While the true overlap area curbs down to a finite value at £2 = 200, (TT2l i diverges. Still, for 
£2 = 201 it produces a fair and for £2 > 202 an excellent approximation. In the degenerate case ( TT3l is indeed capable of reproducing 
the size of the overlap area to good accuracy. 

Thus, we suggest to incorporate degenerate triangle configurations into our formalism by replacing li — > £[ + A^,/2 in all 
arguments of A for these cases. This way, we heuristically correct for the breakdown of approximations in the assignment of the 
actual area, over which triangle configurations are averaged. While the modification is at this stage only motivated by the geometrical 
interpretation, we will establish a more strict foundation of ( TT3l by relating it to Wigner symbols in Sect. 15. II 



4. Bispectrum covariance 

The covariance of the bispectrum is defined as 



The computation of the correlator of two bispectrum estim ators involves a 6-point correlator of g, which can be expanded into its 
connected parts as e.g. outlined in iBernardeau et al. (2002a). Denoting the connected correlators by a subscript c, which will only 
be done in this paragraph to avoid confusion, we obtain 




(13) 



Cov(B(£ u £ 2 ,£3), 5(4,4,4)) = ((M,4,4)-(M,4,4))) (£(4,4,4) - (fi(4,4,4)))) 

= (£(4,4, £3) B(4, £5, 4)) - B(4,4, 4) B(4, 4, 4) • 



(14) 



Wi)# 2 )#3)s(« gits) gih)) 
= (gih) git 2 )) c (gits) gim c (gits) git 6 )) c + (giti) git 2 )) c (g(t 3 ) git 5 )) c (g(t 4 ) git 6 )) c + (13 perm.) 
+ (giti) gih) gih)) c {git a) gits) gih)) c + (giti) git 2 ) g({ 4 )) c (g(£3) gits) git 6 )) c + (8 perm.) 
+ (giti) g(t 2 ) git 3 ) g(t 4 )) c (gits) git 6 )) c + (git^ g(t 2 ) g(t3) g(t 5 )) c (g(t 4 ) g(t 6 )) c + (13 perm.) 
+ (giti) g(t 2 ) git 3 ) git 4 ) g(ts) gih)) c , 



(15) 
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where the permutations are to be taken with respect to the indices of the angular frequencies such that for each correlator, no com- 
bination of indices is repeated (as the individual correlators are invariant under permutations of the indices within that correlator). 
The resulting connected parts are related to spectra via 

N \ ( N \ 

f[ g(jti)\ = (2tt) 2 J] It P N (l { , ... , t N ) , (16) 

i=l ' c \i=l ' 

where we identify P 2 {t\,t 2 ) = P((\) as the power spectrum and P 3 {t\,t 2 , t 3 ) = B((\,t 2 , £ 3 ) as the bispectrum. As the random field 
g vanishes on average, (g(()) = 0, only P4 (the trispectrum) and P^ (the pentaspectrum) will appear in addition in the covariance 
formula, see ( fT5l l. 

Introducing a shorthand notation jf^ ^ d 2 £i/A R ((i) = J., one can write the correlator of the bispectrum estimators by using (f3]i as 

x agto + t 2 + h) ag>(/ 4 + e 5 + 1 6 ) (giA) g{{ 2 ) g (e 3 ) g (t A ) g {i 5 ) g(( 6 )) , 

which then allows us to insert (TT3T > and i\6) . The resulting terms contain products of several delta-distribution. Concerning the terms 
containing three two-point correlators, one obtains e.g. 

5^(t x + e 2 + t 3 ) 5%\t A + e 5 + e 6 ) s^(e x + t 2 ) s%\t 3 + e 4 ) s%\e 5 + e 6 )P(A) P(h) Pits') (18) 

= 8%\e 3 ) 5^(t A +{ 5 + C 6 ) 6^(l x + h) ^{(3 + h) s%\e 5 + P(0) P{€ 5 ) = , 

and likewise for all other terms in which the correlators do not contain one angular frequency each out of the sets {{1,(2,(3} and 
{(4, (5, ((,}■ A similar argument holds for the terms composed of power spectrum and trispectrum, where the trispectrum is readily 
shown to vanish if the two-point correlator contains both angular frequencies out of the same of the sets mentioned above. This way, 
the number of terms with three power spectra reduces to 6, the number of terms with trispectrum and power spectrum to 9. 
To proceed, we demonstrate the treatment of some exemplary terms in the covariance, for instance 

f f f f f f tfftfi +h + (3) ^((4 + ( 5 + 1 6 ) (gih) g (t A )) ( g (e 2 ) g(e 5 )} ( g m 8(h)) (19) 

Jl J2 J3 J4 J5 Jb 

= (inf f f f f f (V D 2) (^ + e 2 + e 3 ) s%\t 4 + e 5 + 4) agto + h) dg>(/ 2 + h) + h) P{t\) Pih) P(.h) 

Jl J2 J3 J4 J5 J6 

A R {(\) A r {£ 2 ) A R (f 3 ) J1J2J3 ' 



where the integrations over £4 to only yield a non-zero result if the annuli of the angular frequencies in the corresponding delta- 
distributions, which are integrated over, coincide. Thus, for every such integration a Kronecker symbol is generated. The resulting 
expression in ( [T9| > can now easily be simplified by using (0, producing a factor of A/(2n) 2 , and subsequently (0 to execute the 
remaining angular integrations. Note that again only the delta-distribution depends on the polar angles of the angular frequencies. 
Therefore, considering only the Gaussian contribution to the covariance, (fTTI i turns into 

(B{€ u iiMB(h,hJ6)) r = - - 7fT~r~7T\ - -j- A" 1 fe,l 2 ,l 3 ) (20) 
\ /Gauss A A R (€x) A R (€ 2 ) A R (€ 3 ) v > 



r Mih C U 2 l 2 C d£ 3 { 3 

x D h h h h h - tb — -J- P(ti) — -J- P{t 2 ) — -J- P(h) , 
Ja s (?j) A r (€i) JA R (e 2 ) A R (i 2 ) JAgfa) A R ({ 3 ) 



where we again pulled A, evaluated at the averaged angular frequencies, out of the radial integrations. Besides, we defined the 
shorthand notation 

Dti&AAAA = ^AA ^<2<5 ^3*6 + hits 8f 2 *t 8*3*6 + 8e,e 4 5e 2 { 6 6( } c 5 + 6( 1 e 5 6f 2 e 6 5(^ t + 6(^ 6 6e 2 t 4 6e } ( 5 + 5c^ 6 6( 2 e 5 6( 3 ( t (21) 
for convenience. By making use of (f2j) and defining the bin-averaged power spectrum as 

M+1/2A« ip f 

Pdi)= _i-ip(^), (22) 

Jl;-\I2M ti&t 

see I Joachimi et al.l ((2008), in analogy to the definition of the bin-averaged bispectrum, one obtains the expression 



(m,€ 2 ,h) £(^5,4)) Gauss - - - - f"J - A"' (*uh,h) DliAAAAA PiWi^Pih) ■ (23) 



A ^^WiA^A^ 
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Terms composed of two three-point correlators can be processed as follows, 



f f f f f f agto + 4 + 4) $o(t4 + 4 + 4) 0K4) g(h) git a)) (g((3) gits) #(4)> (24) 

Jl J2 J3 JA Js J6 

= (2n) A f f f f f fs%\ti + 4 + 4) 8%\t 4 + 4 + 4) <5^ } (4 + 4 + (a) <^ } (4 + 4 + 4) 5(4, 4, 4) B(4, 4, 4) 

Jl J2 J3 JA J5 J6 



= (2tt) 4 6 hh +( 2 + 4) d®(-fi - 4 + 4 + 4) 6%\h + is + 4) 5(4, 4, 14 + 41) 5(4, 4, 4) 

Jl J2 J3 J5 J6 

= A(27r) 2 <^ 4 f f f f ^(4+4+4)^(4+4+4)5(4,4,4)5(4,4,4), 

Jl J2 J3 J5 J6 

where to generate the Kronecker symbol 6[ 1 , we made use of fact that 4 + 4 = _ 4 due to the corresponding delta-distribution. 
To arrive at the last equality, <(5j has been applied after processing the arguments of the delta-distributions similar to ( TT8b . The 
remaining terms, containing four- and six-point correlators of g, can be dealt with in close analogy to (l24b . We mention the special 
case 

[[[[[[ 6 ( o\t 1 + 4 + 4) <4 2) (4 + 4 + 4) 0K4) g(4) g(4)> 0K4) s(4) s(4)> (25) 

Jl J2 J3 J4 J5 J6 

= (2tt) 4 £ f f f f f (^d^ +e 2 + 4)) 2 (^ 2) (4 + 4 + 4)) 2 5(4,4, 4) 5(4, 4, 4) 

= a 2 r r r 4 +4) 5(4,4,4) r r r ^(4+4+4)5(4,4,4) 

Jl J2 J3 J4 J5 J6 

A 2 - - - - 

= A(4,4,4) 4(4,4,4) 5(4,4,4)5(4,4,4), 

which, after inserting this expression into (1171 1. cancels the product B(4,4, 4)5(4, 4, 4) in the definition of the covariance ( fT4l . 
Combining these results, we obtain the total bispectrum covariance 

Cov(b(4,4,4), 6(4,4,4)) = - - - 27r) A _1 f4,4,4) /, a a a. I, 5(4)5(4)P(4) (26) 

v ' A444A4A4A4 1 ; ftAA,WlA 

+ ^ ^ XXXXX^ 2)(/>i + f2 + f3) ^ 2)(/>3 + ^ + #s) B(A ^ 2, B( ^ 3, ^ 4) + (8 perm 
+ ^ ^ XXXXX^ 2 '^ 1 + f2 + f3) ^ 2)(f4 + ^ " ^ 3) p * {tutij * ji) m) + (s perm - 
+ ^ XXXXXX^ 2)(/>i + f2 + f3) ^ 2)( ^ + ^ + 4) p ^ j2j ^ e ^ e ^ e ^ ' 

where the prefactor reads C = (2tt) 6 A" 1 (t x , 4, 4) A -1 (4, 4, 4)- 

The general form of the covariance terms is in agreement with the expressions derived in Sefus attTet alJ {2006). As mentioned 
in Sect.|2] shot or shape noise can readily be included into this covariance by adding a corresponding noise term to the power 
spectra. Weak lensing or galaxy clustering surveys often have in addition tomographic information, so that the data is binned into 
(photometric) redshift bins. The covariance can be generalized to this case in a straig htforward manner by obeying the practical 
rule that each photometric redshift 'sticks' to the angular frequency it is a ssigned to , see lTakada & Jainl d2004l) . A similar argument 
holds for the generalization to CMB polarization bispectrum covariances IHuHTo OO). 

5. Equivalence to spherical harmonics approach 

In this section we demonstrate that both our and the spherical harmonic approach are equivalent in the sense that they measure the 
same information in a survey. Moreover, we investigate the behavior with respect to parity, and the relation between the covariances 
of both approaches, considering for the remainder of this work only the Gaussian part of (1261 1, 

5.1. Comparison of covariances 

On the celestial sphere one can decompose the random field g into spherical harmonics, which produces a set of c oefficients gi m 
with {, m integers and £ > 0, —C < m < t. In terms of the g(,„ one can define a bispectrum estimator as (e.g. Hu 2000) 



) 

) 



ZI i x £ 2 ( 3 \ 
m, m, m, 8e < m > 8e ^ 8t ™ ' (27) 



c u t 2 

m\,m2,mi 



where the object in parentheses is the Wigner-3 j symbol. Properties of the Wigner symbol are reviewed in Hu (2000); most impor- 
tantly, it o beys the triangle condition, i.e. it is non-zero only for 4 - 4 < 4 ^ 4 + 4 an d permutations thereof. For this estimator 
Hu (2000) derived the simple Gaussian covariance 

Cov[Bf u ^ 3 , 5f 4 ,? 5 / 6 ) = -Df,/ 2 ,f 3 ,f 4 ,f 5 ,4 P^P^P^ > (28) 
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where P( denotes the full-sky power spectrum, and where D^um is used as defined in (fJTJ. Moreover, he gives approximate 
relations between the spherical harmonic and Fourier-plane spectra, 



P, - «0 ; B, M . ( << « % ) / M ' + 1)(2 *; lx2f ' + 1 > « , (29, 

valid for t\,lz,£-i » 1. These equations can only hold for integer In addition, the Wigner symbol with m\ — m 2 = m.3 — 
vanishes for L = t\ + ( 2 + £3 odd, see the following section for details. Making use of the standard procedure of multiplying ( f28l ) 
by an ad hoc factor of f^} = 4n/A to account for finite sky coverage of the survey, one can derive a flat-sky spherical harmonic 
covariance with <29) as dHull2000t iTakada & Jainll2004l) 

(B((iJ2, h) B(£ 4 , h, hi) * - (47r) De l M > M ' t l ( \ \ \ \ p(f{) P(i 2 ) P(( 3 ) , (30) 

\ 1 A (2^ + 1) (24 + 1) (2^ + 1) I / 

where still the angular frequencies are required to be integer, and L even. As is true for our approach, (f30b holds for I » 1 only. To 
be able to compare this widely used formula to our results, a relation between the Wigne r symbol and A has to be found. 

When comparing the spherical harmonics and the Fourier-plane approach. iHul d2000h already came across integrals of the form 
(O. We reproduce his computation, 

j d 2 £ { J d 2 £ 2 J d 2 f 3 dg>(/ 1 + h + (3) = J d% j d 2 £ 2 j d 2 £ 3 J ^ e^+^> e (31) 



* J &hU J di 2 { 2 J d^ 3 J|^J" cinY°(n)y£(n)y<?(ii) * 8tt 2 J dA^i j d£ 2 £ 2 j d^ 3 (o q q 

where JdQ is the integral over the unit sphere, and where Y"\ri) denotes the spherical harmonic function with n the unit normal 
vector on the sphere. We are concerned with the validity of this derivation for the following reasons: Terms with integer and real- 
valued I are mixed, e.g. it remains unclear how the integration over the Wigner symbol squared is to be understood. T o get from 
the second to the third equality, the Fourier base e lf is expanded into spherical harmonics, an approximation which IHul ([2000) 
correctly states to be valid for small angles only. However, the integration over angles runs over the full two-dimensional plane or 
the unit sphere, respectively. Moreover, it is not specified how the non-trivial transition from an integral over the plane to one over 
the unit sphere is executed. Instead of OTb . we propose to use (O, which is an exact and rigorous expression. 

To allow for a comparison between (13T) and our approac h based on ©, we need to establish a relation between the square of 
the Wigner symbol and ([§). We refer to lBorodin et al.l d!978l see also references therein) who compute approximation formulae of 
the Wigner symbol in the context of the quasi-continuous limit of quantum states with high angular momenta. The base of their 
derivation is formed by the exact relation 



f dtp f dib f d0smeD e> , (tp, 9, if/) D (l , (tp,6,iJA D h , Up, 9, if/) = 8tt 2 ( £l £l M • ( ^, ^ £ ]\ , (32) 

Jo Jo Jo miin^f' >r/ mim' 2 ^T< <YJ m 3 m, 'V I mj mi mj I \m\ m 2 m\ 

where D , denotes the m x m' element of the Wigner D matrix, which in turn is a function of the three Euler angles tp, 9, and 
if/. Making use of a quasi-classical approximation of the D c mm . ■ iBorodin et aD (1 19781) compute expressions for the general Wigner 
symbol in the limit of large and continuous angular frequencies. From these results we extract the approximation 



l\ t 2 £3 





(33) 



-1/2 



which allows us to generalize the Wigner symbol t o real- valued arguments. Equation d33l holds only for l\, 1%, £3 » 1, which, in 
the quantum-mechanical context of IBorodin et al.l (119781) . originates from the use of expressions that are valid for large angular 
momenta, i.e. the quasi-classical limit, only. This condition on angular frequencies also underlies the approximations in d29b and 
OTI ) and can in our context be interpreted as a natural consequence of working in the flat-sky approximation. 

As is demonstrated in Fig. [4] we find that ( l33l constitutes an exc ellent approximation, whose accuracy over a wide range of 
^-values is orders of magnitude better than the approximation given in lTakada & Jainl (l2004t) . Eq. (A3), 



t\ €2 (3 




,-1/2 ,j ,-1/2 .j ,-1/2 



V2tt 



5-/1 + 1) (--£ 2 + l) (--.3 + 1) (34) 
L/2 -£1 + i/ 2 \ L - 2e,+l/2 (L/2 -£ 2 + ii2\ L - 2h+lil (L/2 -(3 + i/2\^ 2f3+I/2 



L/2 - t x + 1 / \ L/2 - £ 2 + 1 / \ L/2 - £3 + 1 

Only for triangle configurations close to degeneracy does the latter formula perform slightly better. Both approximation formulae 
are least accurate in the case of a degenerate triangle configuration with fractional errors around 10 % or slightly above, but improve 
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10 200 250 300 350 400 10 10000 20000 30000 10 60 80 100 120 140 160 180 

l z ,l 3 triangle area $ 3 [deg] 

Fig. 4. Fractional error of the approximation formulae for the Wigner symbol. Left panel: Shown are the relative deviations of 
( l33l and ( l34l from the true absolute value of the Wigner symbol. The same triangle configurations as in Fig. [2] are used. Results 
for t\ = 200, £3 = 400 and varying I2 are shown in gray while those corresponding to (1,(2 = 200 and varying £3 are plotted in 
black. Solid curves are obtained using (1331 . dashed curves by employing (1341 . Center panel: Same as above, but now plotting on the 
abscissa the corresponding triangle area enclosed by the three angular frequency vectors. Note that before reaching the equilateral 
configuration, the area has a maximum and starts to decrease again. Right panel: Same as above, but now as a function of the internal 
angle <$>?, opposite £3, which is the longest side of the triangle in both configurations considered. Hence, O3 = 60° corresponds to 
the equilateral case, and <J>3 = 180° to the degenerate case. 



quickly to very small percentage deviations when the configuration approaches a more equilateral form. In Fig. [4] we also plot the 
fractional errors as a function of the triangle area enclosed by the three angular frequency vectors and as a function of the internal 
angle <l>3 opposite £3, being the longest side of the triangle in the configurations considered. In terms of these quantities we observe 
a more universal behavior of the errors, in particular in the regime where the approximations are less accurate. We find to good 
approximation that, when approaching the degenerate case, relative errors increase exponentially with decreasing triangle area and 
increasing O3. 

For i » 1, and if the triangle configuration is not too close to the degenerate case, one may approximate {, + 1/2 * C it so that 
one finds from ([8]l and (1331 

t\ h i3f „ A(*i,fc,4l) (35) 



\ ) 2n 

Remarkably, since for integer angular frequencies we have A€\ = M2 = A^? = 1, (1331 exactly reproduces our earlier conjecture 
( fT3] l. which strongly supports its validity. If one replaces the Wigner symbol in OH by d35l ). however, one obtains a result which is 
a factor of 2 larger compared to ©. 

Inserting ( l35l l into (l3"0l . and using 21 + 1 » 2€ for I » 1, we get 

(Wuh, h) HU, h, 4)) * 2n ( ) P(il) P(li) ' (36) 

which is equivalent to (l23l if the latter equation is specified to M\ = A^2 = A^3 = 1, and integer t with L even - except for ( l36l ) 
being a factor of 2 smaller. In the following, we are going to elaborate on this apparent discrepancy. 



5.2. Parity 

To elucidate the different noise properties of the Fourier-plane and spherical harmonic bispectrum estimators, we investigate their 
behavior with respect to parity. In two dimensions the parity transformation corresponds to an axis reflection, or equivalently, the 
reversal of the polar angle of all spatial vectors. To flip the parity of a triangle, one can do an odd permutation of its sides, see e.g. 
the two triangles sketched in Fig.Q] Hence, to test the behavior of estimators for triangles of different parity, it is sufficient to flip 
any two of its angular frequency arguments. 
Consulting d27i >. we find 

K,h,h = (-l/' +f ^ 6 M (37) 
because of the behavior of the Wigner symbol under change of parity, 

i\ ll { } \ _ ,^ti+t 2 +ti I £\ h h J /ggx 

m\ m,2 1113 ) \ mi mi ni2 J ' 
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and likewise for all odd permutations of the columns in the Wigner symbol. Thus, the spherical harmonics estimator is parity- 
invariant for L even and changes sign for L odd. Most cosmological theories predict parity-invariant large-scale structures and CMB 
anisotropies. If parity symmetry is built into the cosmological model at consideration, measures that vary under parity transforma- 
tions do not have any predictive power, wherefore they are usually not considered in a data analysis. Accordingly, ( 127b is only used 
for arguments that have L even. Note that parity invariance is also incorporated into the relation between the spherical harmonics 
and Fourier-plane bispectra, see the second equality of d29l i, via the Wigner symbol which vanishes for L odd (this behavior is a 
direct consequence of d38b for nt\ = mi = 1113 = 0). 

The Fourier-plane estimator is by design parity-invariant, which can be seen mathematically from swapping arguments of (0, or 
illustratively by inspecting Fig.Q] From the sketch it is evident that triangle configurations of different parity are averaged over with 
equal weight. For a more formal argument, we can explicitly construct estimators that averag e only over triangle co nfigurations of 
the same parity. To this end, consider the two-dimensional cross product a x b = a x b y — a y b x (ISchneider & Lombardill2003l) of the 
angular frequency vectors t \, (2, ti- If they form a triangle, one finds t\ X^2 = (2 xtj, = (3 xt 1, which follows from I \ + (2 + (3 = 0. 
A change in the parity of the triangle implies a sign flip in these cross products. 

Noting that t ; x (j = sin(^^ - we compute a condition on the polar angles, 

ipe 2 - (p (i e [0, it] ; tpe, - ipt 2 e [0, it] ; <pt, - (pt 3 e [0, n] . (39) 

To obtain the parity transformed triangle, swap the signs of the polar angles in ( f39b . Under the premise that the vectors do form a 
triangle, one of the conditions in ( 1391 is redundant, the remaining ones restricting the angular integrations in the averaging of ((3). 
For instance, the integration ranges could be modified to <pf t e [0,27r], <pi 2 e [(fe^n + ip^], and (pt 2 € [<p( { - n, (fe, ]. Due to rotational 
symmetry, which still holds, the inner integrals have to yield the same result for all possible values of tpc 1 . Therefore, we can set 
the ranges of the inner integrals to (pt 2 e [0, n] and ip[ 3 e [—it, 0] without loss of generality. To maintain the symmetry, we keep the 
integral over (p^ in our notation. These findings are reflected in the shorthand notation 

d{ip u <p 2 ,ip3] = -r— ( 4 °) 

J Jo 2tt Jo n J_ K n 



which we use to define the following bispectrum estimators, 

f ] -l/2A( J{ 2 -1/2A{ (l^2 J? 3 -l/2Af £3^3 



. 2n 2 u v f ' +I/2Af dl\l\ r^- +ll2M Ahh r 1 ^3(3 

B^e l ,t 1 Ji) = ——A Hi, {2,^3) ■= J ^- (41) 

A J?]-l/2Af tlA£i J( 2 -\/ 2 At €2^(2 Je,-1/2AC 

X d{<f u <f 2 ,<P3] + j d{ip2,?3,<Pl] + j d{(p3,<fi,<P2^j 6 ( u ) ((i+f2 + (3)8((l)g((2)8(e3); 

. 2n 2 u v r { ' +l/2Ae di^ r h - +ll2M d( 2 {2 r h+l ' 

6 v ({ u e 2 ,(3) = — a ((1,(2, £3) ^- L - L 

' Jh-uiM Jh-xnAi ( 2 M 2 Jh-W' 



C,-l/2A{ C1AC1 Jt 2 -\j2Al i-2^2 J{ 3 -1/2A{ 

xijj d{(pu^,n) + J d{<p 2 ,<pu<p 3 } + J d{<p 3 ,<p 2 ,<pi]\ ^{h +t 2 + h)g{t{)g{t2)gm. 

Here, we have symmetrized the restricted integrations d40b by averaging over all either even or odd permutations of [ip^ , ip( 2 ,tpe,}. 
Consequently, changing parity via any odd permutation of the angular frequencies in the arguments of ( f4Tb turns one estimator into 
the other, as demanded, for instance B^(i\,(3, t 2 ) = By{(\,(2, £3)- 

Note that the prefactor of the estimators in d4TT > is diminished by a factor of 2 with respect to (0, which is necessary to keep 
them unbiased. This can be shown by computing the expectation value of (f4TT > in close analogy to the procedure outlined in Sect.ff] 
However, the separate consideration of angular and radial integrals that enabled us to make use of © is not possible anymore in 
this non-symmetric case. For instance, given fixed t\, the restricted angular integrations ( l40t can still produce a triangle of opposite 
parity by including a triangle with \t' 2 \ - {3 and = £ 2 , This is reflected in the fact that the integration (0, if properly normalizecQ, 
still yields the same result when limiting the length of the integration range to n. 

Instead, one can execute the integral over the angular frequency which is still averaged over the full two-dimensional plane, such 



as 



Je l -\!2AC J( 2 -l/2At (2^2 JI 1 -l/2At (3^(3 Jo 27T J 71 J _„ It 

i r c ^ l2hl dhh r e ' +1/2Ae dhh ra^ [° d v± 1 

2itl x hl A }( 2 -n2Ae h^h Jh-i/2U hM 3 Jo n n e ' 



where g f = 1 if (\, ( 2 , (3 form a triangle, and else. The remaining integrations reproduce the overlapping region of the annuli 
for £2 an d {3, as depicted in Fig.Q] By limiting the integration to the half plane to one side of an axis collinear to t\, the overlap is 
obviously halved. Since the area of the annuli for ( 2 and {3 is also reduced by half each, the value of the integration should double, 



1 In the derivation of Sect.|2]the proper normalization of 2n for each angular integral is hidden within A R {1). Note that we have given ((6} without 
this normalization, whereas it is included in i ll It . 
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see dTTT) . Following the geometrical interpretation once again, we thus arrive at 

^22. —12. ^22 _T£L 22ls^\{ 1+ ( 2 +{ 3 )B({i,{ 2 ,{ 3 ) (43) 

Je,-l/2A( €\M\ J( 2 -i/2A£ (2^2 Jf,-l/2Af £3^3 J() 27T J 7T J_ ff 7T 

Comparing this result to (0), the estimators ( l4TT i have indeed to be smaller by a factor of 2 to still be unbiased. 
To obtain bispectrum estimators that are completely analogous to dZTT i. we define 

B ± (li,hJi)= \ {B h {I u hM±Bv{txJ2M) ■ (44) 

As (b-QuhM) 

= for a parity symmetric random field g, and B-{(\,€j,,(2) = -B -(1 1,(2, (3), this estimator shows identical 
behavior compared to B tx with L odd. In practice both measures could be used to assess deviations from parity symmetry. The 
estimators B + {(\,( 2 , (3) and j j with L even are likewise invariant under parity transformations. After some algebra that closely 
follows the outline of Sect.|4]we find that the covariance of B + is the same as (|26| >, which is not unexpected because we already 
noted that (f3]l is also parity-symmetric. 

With (l44l at hand, one can readily extract the different treatment of even and odd parity measures in the spherical harmonic and 
Fourier-plane formalisms. Estimators (|27T i separate the set of possible arguments {(1,(2, (3} disjointly into parity even (L even) and 
parity odd (L odd), whereas B+ and B are defined on the same full set of angular frequency combination^ In other words, when 
limiting B + to integer angular frequencies only, the same information is contained in 'half the number of measures in the spherical 
harmonics case, namely those with L even. The latter estimators have a covariance of half the size of the covariance of B+, so that 
the overall information content is the same for both approaches - as required. 



5.3. Information content 

We veri fy the findings of the foregoing section by comparing the information contained in both approaches in terms of the Fisher 
matrix Qe gmarket al. 1997). For a practical implementation we specialize to a non-tomographic weak lensing survey (see e.g. 
iBartelmann & Schneiderll200ll for an overview), assuming a cosmology-independent covariance that is well approximated by the 
Gaussian approximation, i.e. using d23l and d30l >. respectively. To allow for direct comparison, we limit the Fourier-plane approach 
to integer i with all bin sizes set to unity. Due to the symmetry under permutations of the arguments of the bispectra, one can 
impose the condition i\ < ( 2 < (3 on both formalisms, rendering a block- wise diagonal covariance matrix. Inspecting ( l23l . the only 
dependence on the arguments of the second bispectrum, i.e. (4 to is due to the Kronecker symbols d2Tb . so that the summations 
over £4 to ((, become trivial. 

Hence, the Fisher matrix can be written as 



v dB{( l ,( 1 ,( i ) A (M3 M1M2M3 A(W) dB(( u ( 2 ,(3) 

F>1V = , . . ^aaaa* —af M (2*?mm)m) ' (45) 



where 1 1 » » 1 = 6 for equilateral triangles, D( i j » - t - t » = 2 for isosceles, and Df u t 2 j % = 1 else. The derivatives are 
taken with respect to a set of cosmological parameters p. In this toy example we use only the single parameter Q. m , reducing the 
Fisher matrix to a scalar F. Besides, we restrict the angular frequency values to an unphysically small range between / m i n = 100 and 
!„ = 150 for computational reasons. 

Weak lensing pow er spectra are computed for a standard ACDM cosmology, i ncludin g non- linear evolution via the fit formula 
of ISmith et alj d2003l) . The bispectra are obtained vi a perturbation theory (e.g jFrvlll984T) . using IScoccimarro & Couchmanl (1200 ll) 
with the definition of the non-linear wave vector by Takada & Jain (2004|) to ac count for non-linea r evolution. For the projections 
along the line of sight we assume a redshift probability distribution according to lSmail et al.l ( 1 19941) with B — 1 .5 and a deep survey 
of 0.9 median redshift. Shape noise is incorporated by replacing the power spectra in the covariances with 

P(€) = P{€) + , (46) 
2n 

where the ellipticity dispersion <r e = 0.35 and the galaxy number density n — 40 arcmin 2 are set to typical values for planned 
space-based surveys. 

We calculate the relative deviation of the Fisher information, r = Fp OUI i er /F sp h, harm. - 1, as a function of / max . Note that, since we 
only consider ratios of F, the survey size A drops out. Our results are shown in Fig. [5] For Z max very close to l m [ n one sees alternating 
jumps in r which can mostly be traced back to the fact that, due to the condition l\ < t 2 < €3, the terms entering (l45l l do not always 
split exactly half into L even and odd. After this 'burn in' for / max < 120, r shows only little variation. The remaining offset from 
zero, which is slowly decreasing, can entirely be assigned to the different prefactors in the covariances, i.e. the terms related to 
the Wigner symbol and A, respectively. The range of angular frequencies plotted in Fig.|5]is still far from any physically relevant 
situation, but nonetheless the two approaches agree already better than 99 %. 



2 A similar behavior as for the spherical harmonic estimators would have been unexpected since the possible arguments of B ± form a non- 
countable set. 
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Fig. 5. Comparison of the Fisher informa- 
tion as obtained by spherical harmonics and 
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6. Conclusions 

In this work we intended to give insight into the derivation and the form of the bispectrum covariance in the flat-sky approximation, 
based exclusively on the two-dimensional Fourier formalism. We defined an unbiased estimator that takes the average over the 
overlap of annuli in Fourier space, and computed its covariance. To obtain precise normalizations, a case distinction is necessary 
between degenerate and non-degenerate triangle configurations. However, given that both normalizations become very similar for 
{ » 1, which is assumed in the flat-sky approach anyway, we suggest as a simple and fair approximation to use the expression 
derived for the degenerate case. Then our result for the Gaussian part of the bispectrum covariance reads 

Cow(B{£ u e 2 ,h), B(l A ,e 5 ,lS\ = (2 - 7r _ ^MlMlMl A" 1 [l x + —,€ 2 + —J 3 + —) Wi)P(4W 3 ) . (47) 

This formula is readily generalized to the total covariance by modifying the arguments of A, appearing in the non-Gaussian terms of 
0261 1. accordingly. It is directly applicable to any real values of angular frequencies, to arbitrary binning, and to any compact, finite 
survey geometry. This formula can be modified to incorporate shot or shape noise, as well as to account for photometric redshift 

information or CMB polarization in a straightforward manner. 

While the general form of our result was in agreement with existing work, we found, contrary to I Iu (2000), that the size 
of the covariance is a factor of 2 larger than the one obtained by the flat-sky spherical harmonic approach. By defining parity- 
sensitive bispectrum estimators, we discussed the behavior of both formalisms with respect to parity transformations, arguing that 
the difference in the covariances is indeed to be expected because in the spherical harmonic framework, parity-invariant measures 
are restricted to a subset of the angular frequency combinations at which the bispectra are evaluated. In a practical example we 
demonstrated that both approaches indeed contain the same information in terms of the Fisher matrix, with a h igh level of agreemen t. 
As a consequence, we can confirm that studies performed in the flat-sky spherical harmonic approach, such as lTakada & Jainl (|2004), 
yield correct parameter constraints as long as the analysis is restricted to integer I with the sum of the three angular frequencies 
being even. 

We established a relation between the geometrical and intuitive process of averaging over the overlapping regions of annuli in the 
Fourier plane and the Wigner symbol of the spherical harmonic approach. Both quantities were demonstrated to be in turn connected 
to a simple measure that is proportional to the size of the area enclosed by the triangle configuration for which the bispectrum is 
calculated. This resulted in convenient, yet precise approximation formulae for the prefactors of the covariances of both approaches 
at consideration. 

Under the assumption of a compact survey geometry and scales much smaller than the extent of the survey area, (l26l i provides 
a cleanly derived bispectrum covariance matrix that naturally incorporates the scaling with survey size, is not restricted to integer 
angular frequencies, and allows for any appropriate binning. 
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